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The dynamic behavior of cluster algorithms is analyzed in the classical mean held limit. Rig- 
orous analytical results below T c establish that the dynamic exponent has the value z aw = 1 for 
the Swendsen - Wang algorithm and z U w = for the Wolff algorithm. An efficient Monte Carlo 
implementation is introduced, adapted for using these algorithms for fully connected graphs. Ex- 
tensive simulations both above and below T c demonstrate scaling and evaluate the finite-size scaling 
function by means of a rather impressive collapse of the data. 



I. INTRODUCTION 

In a series of seminal pioneering papers Fortuin and 
Kasteleyn |ij established an important connection be- 
tween percolation problems and the related concept of 
clusters on the one hand and various thermodynamic 
functions of Potts spin models on the other. These con- 
nections allow geometric interpretation of different prop- 
erties of these systems, that play a central role in their 
critical behavior. For example, the onset of spontaneous 
magnetization coincides with the appearance of an infi- 
nite cluster in the related percolation problem; the sus- 
ceptibility is proportional to the mean cluster size, and 
the correlation length to a typical cluster's radius. 

Some years later Swendsen and Wang (SW) Q used 
these ideas in a way that constituted a breakthrough in a 
different sub-field of statistical physics: that of computer 
simulations of systems in thermal equilibrium. They ex- 
ploited the Fourtuin-Kasteleyn mapping to define a very 
efficient Monte Carlo algorithm that performs large scale 
non-local moves, by flipping simultaneously entire clus- 
ters of spins. This dynamics decreases the relaxation 
time to equilibrium while preserving detailed balance and 
ergodicity. In this way one can overcome (or at least 
significantly reduce) critical slowing down at second or- 
der phase transitions ||. "Critical slowing down" is a 
well known phenomenon; the relaxation time of a sys- 
tem diverges as the critical point is approached. The 
manifestation of this real, physical effect on simulations 
is that at criticality the typical time needed to produce a 
large set of decorrelated configurations, that appear with 
the proper Boltzmann weight, diverges with the system's 
size. Hence the standard local Metropolis-type Monte- 
Carlo simulation methods become inefficient. 

A wide variety of cluster methods have been applied 
successfully to many fields in physics (second-order phase 
transitions, disordered and frustrated systems, quantum 
field theories, fermions in a gauge background, quantum 



gravity and more p|-pl|). In many cases a dramatic accel- 
eration of convergence to equilibrium was achieved. That 
is, if the simulation of a system of (linear) size L is per- 
formed at the critical point, the relaxation time behaves 
according to 

where z c , the exponent associated with the cluster algo- 
rithm is significantly less than zi oca i of the local methods. 

So far only relatively few rigorous results have been 
derived for cluster algorithms and the associated relax- 
ation times and exponents |^,^| . In this paper we consider 
the well known version of the Ising spin model for which 
mean field is the exact solution; that of N — ► oo spins, all 
connected to each other. By analyzing this system one 
can calculate the classical, mean-field value of z c . 

The mean-field Ising ferromagnet has been previously 
examined analytically and numerically by Ray, Tamayo 
and Klein (RTK) 0J. Although the work presented here 
agrees with the conclusions of RTK regarding the SW 
dynamics, it contains some improvement. Firstly, the 
derivation of the analytical results is improved, employ- 
ing the known exact results of Erdos et. al. for 
highly diluted random graphs, instead of the approxima- 
tion (used by RTK) of a Bethe lattice with large coordina- 
tion number. We show that this technical improvement 
does not alter the results obtained Q for SW dynamics. 
Secondly, we performed simulations, introducing a new 
efficient algorithm, which extend those of RTK in the 
following aspects. RTK looked at the relaxation times r 
of the magnetization autocorrelation function, measured 
at a single temperature (T c of the infinite system) versus 
the size of the system. In addition they also studied the 
temperature dependence of r below T c at a single value 
of the system size. Our simulations were done on larger 
systems using longer running times; in addition, we stud- 
ied a wide range of system sizes and temperatures both 
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above and below T c . Our analysis of the results is based 
on finite size scaling; the relaxation times measured at 
various temperatures and system sizes are presented us- 
ing data collapse, thereby also calculating numerically 
the corresponding finite size scaling function. This data 
collapse of course takes into account the shift in the tran- 
sition temperature as a function of the size of the system. 
It resolves the unexpected result which RTK found puz- 
zling, namely that the slowing down below the transition 
appeared to be enhanced in comparison to that at the 
transition. When the shift in the transition temperature 
is properly accounted for, this effect disappears. 

In addition, we present analytical and numerical re- 
sults for Wolff (UW) (n[] dynamics with the conclusion 
that the relaxation times of the autocorrelation function 
depends on the details of the global algorithm. However, 
the upper critical dimension is hypothesized to be unaf- 
fected by the details of the local or the global algorithm. 

II. THE MODEL 

The Hamiltonian of the fully connected Ising model is 
defined by 

1 N 

The statics of this model is solved exactly, using mean- 
field [l3| |. The result is a second order phase transition 
at p c = 1, with exponents a — 0, v = 0.5, (3 = 0.5 and 
7 = 1. The local (single spin flip) dynamics of this model 
can also be solved exactly using a mean-field approach re- 
sulting in (from here we use (3 = 1/kT unless otherwise 
stated) 

t(13)k((3-[3 c )-\ 

Using the fact that v — 0.5 (i.e. £ oc (/3 — (3 C )~^), it is 
clear that toc( 2 , giving 

Zlocal — 2. 

The standard SW procedure as implemented for this 
model is defined as follows. Let us assume that the sys- 
tem is in a spin configuration {Si}, i — 1..N. In each SW 
sweep every bond is "activated" (or frozen) with proba- 
bility 

P f = l-exp{-^(l + S i S j )}. (2) 

A bond that is not frozen is called "deleted". Frozen 
bonds define connected clusters, to each of which a spin 
value (±1) is now assigned randomly, giving rise to the 
new spin configuration. Then the procedure is repeated. 
The dynamics of this SW procedure was first studied for 
the fully-connected model by Ray et al who showed 
that z sw = 1. 



We present now a different derivation of this result - 
one which utilizes exact results from graph theory. 

In order to get better insight into the percolation prop- 
erties of the fully connected model let us introduce a few 
previously known facts about highly diluted graphs [fl2f . 
These are relevant to our problem since for N >> 1 the 
freezing probability Pf is very small, and the SW pro- 
cedure generates precisely the highly diluted graphs to 
which we now turn. 

Consider a graph of N nodes; on all links connecting 
pairs of nodes i,j place independent random variables 
Jij, taken from the probability distribution 

P(J t3 ) = (1 - c/N)5{J tj ) + (c/N)6(J ZJ - 1) (3) 

where the average connectivity c is taken to be 0(V). 

The resulting highly diluted graph has finite connec- 
tivity; in the thermodynamic limit the probability that 
a node has connectivity k (i.e. has non-vanishing link 
variables Jy with k other nodes) follows the Poisson dis- 
tribution: 

P(k) = c fc exp(-c)/fc!. 

Many geometr ic p roperties of this diluted graph are 
well understood |]l4p^ , ^5| . In particular, the graph un- 
dergoes a percolation transition at c = 1. That is, the 
number of nodes that belong to the largest connected 
cluster is of order (log TV) for c < 1, 0{N 2 / 3 ) at c = 1, 
and of 0(N) for c > 1. In this percolating regime the 
size of the largest cluster is gN, where the parameter g 
is the solution of the equation |Q 

.9=1- e- cs 

or 

ln(l - g) = eg. 

Assuming that for c = 1 + 6 and 6 <C 1 one has g <C 1 we 
get 

cg-g- g 2 /2 - g 3 /3 = 

and therefore 

2g 2 + 3g- 66 = 

which yields 

g(l + 6)=26-^. (4) 

After this aside on random dilute graphs we return 
to the SW dynamics for our fully connected Ising spin 
model. Since the spins are fully connected, a configura- 
tion is fully parametrized by N + , the number of positive 
valued spins, or by its magnetization 

m = (N+ -N-)/N 
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where N_ = N — N+. The dynamics is, therefore, com- 
pletely characterized by specifying the value of either iV + 
or m after each Monte Carlo sweep. Every SW sweep 
can be thought of as consisting of two steps. In the first 
step the spins are divided into two groups, one of positive 
spins (all the spins with Si = + 1) and the other negative 
(spins with Si = —1). Using Pj of (g) we first note that 
all the links connecting any positive spin to any negative 
one are deleted. The positive (negative) group contains 

N 

N ± = -(l±m) 
spins. Without loss of generality, we assume that N + > 

In the second step, each group is examined indepen- 
dently. The internal links within a group are " activated" 
with the probability Pf = 1 — exp(—2(3/N). Each group 
can be thought of now as a random graph with effective 
connectivity 



c± = P f x N± 



2W± 
N 



(5) 



Denoting by t the deviation from criticality: /3 = 1 — t 
(recall that (3 C = 1) we can write for the larger (positive) 
group of spins 



c+ = (l-t)(l + m). 



(6) 



We can now view the act of freezing a bond as choos- 
ing Jij = 1 in the generation of the dilute graph dis- 
cussed above. For that problem, as was mentioned pre- 
viously, c = 1 is a critical value. Therefore in the disor- 
dered phase of our spin model t > and j3 < (3 C so that 
both c_|_, c_ < 1 and the SW freeze-delete procedure cre- 
ates only small (non-macroscopic) clusters. On the other 
hand, at temperatures just below criticality, i.e. when 
t < 0, we have c+ > 1 but c_ < 1. Hence only one 
macroscopic cluster (of positive spins), whose size is pro- 
portional to the lattice size, is generated, while the rest 
of the clusters are small. 

Flipping of this single macroscopic cluster is the domi- 
nant element of the dynamics. The rest of the clusters are 
small, there are 0(N) of them, and therefore the extent 
that their flip influences the magnetization is negligible 
compared to that of the big cluster. 

In what follows we derive recursion for N l + and mi, 
the values that the number of the majority spins and the 
magnetization take after I steps of the SW procedure. 
We use the notation N + , but mean the number of spins 
in the state with majority (even when they are negative). 

The small clusters become ±1 with probability 1/2 af- 
ter a sweep: neglecting fluctuations of their contribution, 
and assuming that the single large spanning cluster was 
and remains positive, we can write 

< +1 = \Nl + \[N l + - N l + g(c + )} + N l + g(c + ) (7) 



N l _ +1 = -N l _ + -[N l + -N l + g(c + )}. 



(8) 



Denoting ni = N l + /N equations (7-8) can be rewritten 
as a recursion for the magnetization: 

mi+i = n ig {2f3n[) = -±^<?[(1 - t)(l + mi)}. (9) 

For simplicity we keep m positive (or, in other words, our 
recursion is actually for |m|). 

Near the transition the argument of the function g is 
near its critical value and we can identify <5 (see Eq. (4)) 
by 

5 = — t + mi — tmi. 
Equation (4) is now given (up to terms of order t 3 / 2 ) 



by 



1 + mi , 8/2 m 

mi+i = — ^ — i 2 v * + mi ~ tmi > ~ 3V 771 ; _ tm i)\ 



2 1 2 

= mi - t + -mit - -m l , 

which leads to the differential form: 



dm 



-t H — md 
3 ' 



(10) 



(11) 



where I, the discrete index of the sweep has become a 
continuous "time". For t < Eq. ( pi] ) has a stationary 
point at 



m* = (-34)3. 

Near this point we write 

m = m* + e, 

and linearize equation ( |Tl| ) in e: 

de 2 2 . 2 2 . 

— = -et m,e = -et ■=e(—ty- 

dl 3 3 1 3 



(12) 



(13) 



For \t\ <C 1 the (— t)? term is dominant so one can con- 
clude that relaxation is dominated by 



e ex exp 



and therefore in the ordered phase we find the relaxation 
time 



>/3, 



(14) 



Using the definition 



r =(-*)" 



yields our result for the mean field value of the dynamic 
exponent of the Swendsen-Wang procedure: 

•£.s?(? I- 
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III. NUMERICAL RESULTS 
A. The Algorithm 

The standard version Q of the SW algorithm vis- 
its, during a single sweep, each bond once and takes a 
freeze/delete decision 

The fully connected system (Eq. 0) contains N spins 
and N 2 /2 links; consequently the usual SW algorithm 
requires 0(N 2 ) operations for one sweep of the system. 
For our choice of Pf (Eq. |2|), however, only a very small 
fraction of bonds is eventually frozen: even though every 
spin has 0(N) links, only 0(1) are activated. Moreover, 
as explained in the Introduction, in the fully connected 
model the only relevant parameter that determines the 
system's state is its magnetization, or N + . Hence to get 
the new configuration after the SW sweep, all we have 
to determine is the new value of N + . To do this we do 
not have to scan all O(N) bonds of a certain spin; rather, 
when the algorithm "visits" a new spin that belongs to a 
cluster, it suffice to use the binomial distribution to de- 
termine the number of activated bonds it has (i.e. the 
number of new members of the cluster). 

Therefore the SW algorithm can be implemented for 
the fully connected model by manipulating a few num- 
bers, instead of the values taken by the spins. This 
approach simplifies the algorithm, is much more effi- 
cient, and uses a very small memory. 

We adopted the spirit of the Wolff algorithm to grow 
all SW clusters during a sweep. At a given moment we 
know the number of spins that have been assigned to 
our presently growing cluster. New spins that have been 
added to the cluster have to be tested for the number of 
frozen bonds connecting them to free spins (that do not 
yet belong to any cluster and have the same sign as the 
presently growing cluster and hence are candidates for 
joining it). The spins that are known to belong to the 
cluster but have not yet been tested are in a "stack". As 
a cluster is growing, the algorithm tracks and modifies 
the following parameters: 

• N c = number of spins in the growing cluster 

• N s = number of spin in the stack 

• Nf — number of "free spin" 

A SW sweep starts at a configuration with N + " up" spins 
and N_ " down" . Say we start to generate clusters of the 
up spins. The following steps are taken: 

1. Assign one spin to the growing cluster: set N c — 
N s = 1, Nf =N+-l. 

2. Test the spin: that is, using the binomial distribu- 
tion, determine N add , the number of frozen links 
our spin will get, connecting it to presently free 
spins: 



(15) 

3. Update : 

• N c = N c + N add , 

• N s = N s + N add - 1, 

• N f = N f - N add . 

4. Perform another test (as long as N s > 0) using step 
#2. Repeat 2 — 4 until the stack is empty (N s = 0). 
construction of the cluster has now been completed. 

5. Flip the cluster with probability \. If the cluster 
has flipped, the value of the magnetization changes: 

M — > M — 2N C . 

6. Reset N c = 1, N s = 1, Nf = N f - 1 and repeat 
2 — 6 to construct new clusters, until there are no 
more free spins (Nf = 0). 

7. Go to the second component (of down spins) and 
apply 1 — 6 with N-. Note that now when the 
cluster flips the magnetization increases by 2N C . 

Note that all the operations of our algorithm arc on a 
few numbers, namely we do not really define spins, just 
manipulate a few counters that were defined above. 

For the Wolff version one can apply the same algorithm 
and grow only the first cluster; the decision of starting 
it from an up or down spin as seed is done at random 
(with probability proportional to the number of up/down 
spins). 



B. Finite Size Scaling Analysis 

Finite size scaling is a standard method |[7|,[l8| to ex- 
tract critical indices from numerical data obtained by 
simulating finite systems of linear size L. Denote by t 
the reduced temperature variable, and by B a thermo- 
dynamic quantity that exhibits singular behavior in the 
L — > oo limit, characterized by the exponent 9: 

Boer 9 . (16) 

For a finite system clearly B is an analytic function of t 
and the finite size scaling function exhibits a crossover; 
when t is such that the correlation length £ ~ t~ v is 
much smaller than the system size L the function be- 
haves like (^6[); in the L — > oo limit all dependence on 
L should disappear. On the other hand, for £ < L the 
finite size effects wipe the singularity out. 
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The scaling hypothesis [jl7| asserts that the crossover 
scaling function depends on L only through the ratio L/£ 
and one expects the form 



B< k t 1 L) = L B ' v f{LIZ) = e' v g{LlZ). 



(17) 



This scaling form includes implicitly the finite size shift 
of T C (L) which defined by the temperature at which the 
measured quantity has a peak. The size of this shift gen- 
erally has the following form: 



|T c (i) - T c (oo)| ocL 1 ^. 



(18) 



The specific scaling functions which we consider in this 
paper are: 

For the susceptibility: 



and for the relaxation time: 



-(t,L) = L'f(tL 1 ' v ). 



(19) 



(20) 
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FIG. 1. The susceptibility, \, versus (3 for different system 
sizes, measured in the vicinity of the bulk critical temperature 

I3 C = 1. 

The most striking manifestation of scaling is data col- 
lapse. The product \li~il v should depend only on 
L/£ f» A/3, where A/3 = — 1. This can be seen 
very clearly in Fig. 2, where the mean field exponents 
7/^ = 2,^ = 1/2 were used (we plotted the scaling 
function versus (L/£) 2 ). 



C. Details of the Simulations 



The Monte Carlo simulations were carried out with 
the algorithm described above, for systems that contain 
N = 10 4 — 3 x 10 5 spins. 

The number of MC sweeps per spin was 10 5 . The first 
10 4 configurations of the simulation (at list 500 times the 
relaxation time itself) were discarded in order to ensure 
that the system reached equilibrium. 

The static equilibrium property that we measured was 

X = f3N(< m 2 > - < \m\ > 2 ), 

which is proportional to the "real" susceptibility x 0- 
The reason for this choice is that in finite systems, espe- 
cially with cluster dynamics, the quantity < m 2 > — < 
in > 2 has no peak at all and the measured value in- 
creases monotonically as T decreases. The reason is that 
m changes sign frequently and therefore < m >^ 0. The 
method that determines x by measuring the size of the 
SW clusters jnj also yields no peak. 

Fig 1. shows the susceptibility for different system sizes 
L as a function of temperature. One can see the large 
shift of the peak position for different sizes of the system. 
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FIG. 2. The susceptibility scaling function (Eq. (|19|)). The 
mean field value of j/v = 2 was used for collapsing the data 
of Fig. 1. 



Besides presenting the scaling function and verification 
of the scaling form, the excellent data collapse that can 
be seen verifies the accuracy of our numerical simulations. 

The shift of the effective T C (L) is studied in Fig. 3; 
indeed one can see that the size of the shift agrees very 
well with Eq. [l^, with the theoretical mean field value 
v=\. 
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FIG. 3. The shift of the peak Af3 max = (L) — 1 of 
the susceptibility x> an d °f the integrated relaxation time r, 
plotted versus L — N%. Non- linear least square fit yielded 
the values l/v = 1.98 ± 0.12 for x, and 1/v = 2.01 ± 0.05 for 

T. 



The dynamic property that we measured is the time- 
dependent correlation function of the absolute value of 
the magnetization. The measurements were taken after 
every sweep. We considered the normalized correlation 
function of the absolute value of the magnetization 



Pit) 



< \m(t')m(t' + t)\ > - < \m\ > 2 > 
< m 2 > - < Iml > 2 > 



(21) 
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FIG. 4. The integrated relaxation time t, plotted versus (3 
for different system sizes, in the vicinity of the bulk critical 
temperature f3 c = 1. 

For resonable sizes of the systems (such as the sizes 
that could be treated with computational power avail- 
able to us) one can not neglect this shift. For example, 
calculating r (as well as any other quantity) for different 
system sizes at T c (oo)(= 1) and trying a fit of the form 

r(T c (oo))=i*, 

introduces some systematic errors. 

Fig. 5 presents the first measurement (to our knowl- 
edge) of a scaling function for r. The theoretically de- 
rived value z = 1 was used. 



For which the integrated relaxation time, defined by: 



Tint 



-. oo 



(22) 



was calculated. 

Fig. 4 shows r for different sizes of the system at vari- 
ous temperatures. As in the case of the susceptibility, one 
can see the large shift of the peak of the scaling function 
for different system sizes. Fig. 3 contains also the scaling 
analysis of this shift as a function of L: again, the result 
agrees very well with the theoretical prediction. Note 
that although the slope is the same as in the suscepti- 
bility case the shift is much bigger, i.e. the amplitude is 
larger. 
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FIG. 5. The scaling function for the integrated relaxation 
time (Eq. ©). The value z = 1 was used for collapsing the 
data of Fig. 4. 

The evident data collapse indicates a good fit to the 
scaling form presented in Eq. (^0|) and to the predicted 
z. 
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IV. WOLFF-DYNAMICS 

Let us start with a remark about SW dynamics. Going 
back to Eq. || one can see that for a particular version 
of SW algorithm, in which one flips only the spanning 
cluster, the following equations are obtained 



N l+i 



N i+i = N i_ + N l + g{c+). 
Denoting N l + /N by ni, this yields 
m;+i = mi - 2nig(2(3ni) 



(23) 



(24) 



mi — (1 + mi)g 



(1-0(1 + 771;) 



(25) 



Keeping m positive, one would get instead of Eq. [Tl| 



dm 



-2t 



-mit 



(26) 



Going over the previous calculation (Eqs. [Tl]-p~4t) one 
gets the same fixed point, the critical exponents do not 
change, and only a factor ^ multiplies the prefactor in 
the expression for r. That is, 



This result shows that flipping the spanning cluster con- 
stitutes the central component the relaxation process. 
Moreover, it shows that flipping small clusters together 
with the flip of the spanning cluster does not help the 
relaxation at all! 

We turn now to UW dynamics and calculate t uw (/3) in 
the thermodynamic limit, for T > T c . We calculate t uw 
in two steps. In the first step we calculate r"™, the re- 
laxation time measured using every single cluster update 
as one unit of time. In the second step, we normalize 
these UW time units to regular time units (N spins are 
accessed or updated in a single regular unit of time) . To 
achieve this normalization we multiply the time measured 
in UW time units by the ratio of the average cluster size 
< C > to the system size N. 

To calculate r! 1 ™ note that 



N 



c 
1 

: - x — 
9 N. 



1 . 
-r 
2 

N 



N 

Z 7 



(1 + m)g 



—, (27) 



9 



where C — gN is the size of the macroscopic cluster. To 
understand the leftmost equality note that every cluster 
update of the UW procedure has a chance of % to hit 
the macroscopic cluster. The rest of the cluster updates 
are irrelevant in the thermodynamics limit. When nor- 
malized to regular units of whole lattice update we get 



< C > 
N 



(28) 



where < C > is the average cluster size in the UW pro- 
cedure. 

We turn now to calculate < C >, using the theory of 
random graphs jjq ]. For the part c < 1, and the clus- 
ter size distribution has only "small-clusters" with aver- 
age cluster size 7*1 =0(1). For the iV + part c > 1 and the 
cluster size distribution has two components. The first 
component is that of " small-clusters" with average clus- 
ter size r 2 = 0(1) and the second component is a single 
spanning cluster, whose size depends linearly upon JV+ 
with size g{c+)N + . Therefore: 



< C > 
N 



1 

N 



Iv 



N. 



■TV 



N 1 ^ 

w— cluster size£N+ 



P(w) x w} 



1 

N 



+ ^±{(1 - g{c + ))r 2 + g(c + )(g(c + )N + )} 



In the limit N —* oo, the first two terms vanish, and the 
leading term is (~]^r) x 9' ''■ Thus, we arrive at : 



< C > _ (1 + mfg 1 ^ g 2 



N 4 
Combining equations (|b4]),(p7|) and (p8|), we find that 



(29) 



t g(c+){l + mJ/4 = r g 



(l-t)(l + m) 



(l+m)/4, 
(30) 



Using equation (||) we expand g to obtain: 

T i 

V3, 



T uw sa T sw 2(-t + m - tm)(l + m)/4 sa r sw m/2 



'-(-t)-i x V3(-t)i/2 w 0.75 = Const. (31) 

Thus we have shown that t uw does not depend on t and 
hence z = for T < T r . 



V. NUMERICAL RESULTS- WOLFF 

The Monte Carlo simulations were carried out using 
the algorithm described in Sec 3, growing only the first 
cluster. Systems that contain N — 2 x 10 4 — 3 x 10 5 
spins were simulated. The number of cluster flips was 
2 x 10 6 for each run. The first 10 5 configurations of each 
simulation were discarded in order to ensure equilibra- 
tion. As in the SW case, the dynamical property mea- 
sured was the correlation function of the absolute value of 
the time-series of the magnetization. The measurements 
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were done after every 4 cluster updates. ^From the cor- 
relation function we calculated the integrated relaxation 
time. Finally we normalized by the average cluster size 
as was measured during the simulation. 

Fig. 6 shows t for different sizes of the system at vari- 
ous temperatures. In order to get the scaling function we 
performed data collapse for r with the theoretical value 
of z = 0. 
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FIG. 6. The Wolff integrated relaxation time t uw , plotted 
versus /3 for different system sizes, in the vicinity of the bulk 
critical temperature /3 C — 1. 



VI. SUMMARY AND CONCLUSIONS 



Cluster algorithms constitute a significant advance in 
computational physics. Their associated dynamics is 
characterized by critical exponents that differ from those 
of standard methods. Very few exact results have so far 
been derived regarding these exponents. 

In the present work we addressed the issue of the clas- 
sical limit for the dynamic exponents of two cluster al- 
gorithms. For the Swendsen-Wang algorithm we found 
z = 1; our analytical derivation of this result is based 
on rigorous results from the theory of random graphs. 
Since the proof applies only for the low-temperature 
phase, we measured the relaxation time t sw , using a 
novel implementation of the SW algorithm adapted for 
fully-connected graphs, both above and below T c . We es- 
tablished scaling behavior of t sw (t, L) by demonstrating 
rather conclusive data collapse. Previous work used a 
Bethe-lattice approximation for the analytic derivation, 
and did not present the full scaling function for the re- 
laxation time. 

For the single-cluster algorithm of Wolff we have shown 
that z = 0, and presented numerical evidence, again ob- 
tained both above and below T c . The scaling function 
was found numerically and data collapse was demon- 
strated. 
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The results is shown in Fig. 7, in which one can see a 
very good agreement with the scaling form presented at 
Eq. In addition the asymptotic behavior is consistent 
with the result of Eq. ([H]), valid for (3 > (3 C and L » I 
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